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CONTINUITY OF A QUEUEING INTEGRAL REPRESENTATION 
IN THE Ml TOPOLOGY^ 

By Guodong Pang and Ward Whitt 

Columbia University 

We establish continuity of the integral representation y{t) = x{t) + 
h(y{s)) ds, t > 0, mapping a function x into a function y when 
the underlying function space D is endowed with the Skorohod Mi 
topology. We apply this integral representation with the continuous 
mapping theorem to establish heavy-trafSc stochastic-process limits 
for many-server queueing models when the limit process has jumps 
unmatched in the converging processes as can occur with bursty ar- 
rival processes or service interruptions. The proof of A/i -continuity is 
based on a new characterization of the Mi convergence, in which the 
time portions of the parametric representations are absolutely con- 
tinuous with respect to Lebesgue measure, and the derivatives are 
uniformly bounded and converge in Li. 

1. Introduction. The integral representation 

(1.1) y(t)=x{t)+ [ h{y{s))ds, t>0, 

Jo 

mapping a function x into a function y, plays an important role in heavy- 
traffic stochastic-process limits for many-server queues. Theorem 4.1 of our 
review paper [4] shows that this integral representation constitutes a con- 
tinuous map on the function space D = D([0,T],M) with either the uniform 
or Skorohod Ji topology [7], provided that the function /i : M — t- M appearing 
in the integrand is a Lipschitz function, that is, 

(1.2) \h{wi) — h{w2)\ < c\wi — W2\ for all wi, tt;2 € M, 

with c being the Lipschitz constant. As a consequence, the integral repre- 
sentation can be applied with the continuous mapping theorem to establish 
desired stochastic-process limits. 
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Our purpose here is to extend that continuity result to the Skorohod Mi 
topology [7], as discussed in Chapter 12 of [8]. As illustrated here in Section 
2, that enables us to obtain associated stochastic-process limits when the 
limit process has jumps unmatched in the converging processes (see Chapters 
1 and 6 of [8] for additional discussion). The desired result is the following 
theorem. 

Theorem 1.1 (Continuity in Mi). The function ip-.D^D mapping x 
into y, defined by the integral representation (1.1) with h Lipschitz as in 
(1-2), is continuous if D is endowed with the Mi topology. 

In order to establish this result, we develop a new characterization of 
convergence in (Z),Mi). Since this characterization is likely to have other 
applications, it is of interest in its own right. Indeed, the bulk of the paper 
is devoted to this extension. To state the result, recall that x„ — )■ x in D if 
dMi {xn,x) — )■ as ?i — )■ oo where dj/j is the metric 

(1.3) dMAxn,x)= ^ inf - n|| V ||r„ - r||}, 

V if;2 = max{t(;i, ^2} for wi,W2^^, = supo<^<i{|u(s)|} and n(x) is 
the set of all parametric representations (u,r) of x. 

A parametric representation [u, r) is a continuous nondecreasing function 
of the interval [0, 1] onto the completed graph of x, where the function 
u gives the spatial component, while the function r gives the time compo- 
nent. In this context, "completed" means that the graph contains the sets 
{{px{t) + (1 — p)x{t—),t) :0 <p < 1}, so that the graph is a connected com- 
pact subset of M X [0, T] , while "nondecreasing" is with respect to the order 
following the continuous path on the graph in M? starting at (with infimum) 
(a;(0),0) (see [8], page 81, for more details). 

As indicated by Theorem 12.5.1 of [8], there is considerable freedom in 
the choice of parametric representations. We will want to use versions such 
that the time components r„ and r are absolutely continuous with respect to 
Lebesgue measure and have uniformly bounded derivatives, where there is Li 
convergence of the derivatives as well as convergence of the time components 
themselves, as in (1.3). For that purpose, let the Li norm of the function r 
be 

||r||i = / \r(s)\ ds. 



Theorem 1.2 (Time functions in the parametric representations). Sup- 
pose that — )• X in {D,Mi) as n— t-oo. Then there exist {u,r) and {un,rn) 
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as parametric representations of x and Xn, where both r and r„, are abso- 
lutely continuous with respect to Lebesgue measure on [0, 1] with derivatives 
r' and r'„ for all n such that 

(1.4) 11?'^, — ''''111 — ^ as n — 7- cx), ||r'|| < oo and sup{||r[j||} < cx). 

n>l 

Our proof of Theorem 1.1 is based on a simple change of variables, much 
like the Ji argument in [4]. For the Mi topology, it exploits the structure 
provided by Theorem 1.2. 

The rest of this paper is organized as follows. In Section 2 we apply 
Theorem 1.1 to establish a many-server heavy-traffic stochastic-process limit 
for the G/M/n + M model where the scaled arrival process converges to a 
limit with jumps but with modified scaling. In Section 3 we prove Theorem 
1.1; in Section 4 we prove Theorem 1.2. 

2. Many-server heavy-traffic limits with unmatched jumps. We now ap- 
ply Theorem 1.1 to obtain heavy-traffic stochastic-process limits for many- 
server queues. Our result here extends previous limits for the G/M/n + M 
model in Theorem 7.1 and Section 7.3 of [4] to cover nonstandard scaling and 
jumps in the limit process that are unmatched in the converging processes. 
For an application to treat service interruptions in many-server queues, par- 
alleling earlier work by Kella and Whitt [3] for single-server queues; see [5]. 

Here we consider a sequence of G/M/n-\-M queueing models with general 
arrival processes (the G) and customer abandonment (the +M), indexed 
by the number of servers, n. For each n > 1, the n homogeneous servers 
have independent exponential service times with rate and customers have 
independent exponential patience times with rate 6. 

Let the arrival rate in model n be and assume that Xn/n — t- A > as 
n — >■ oo. Let An{t) count the number of arrivals in the interval [0,i]. We 
assume that the arrival processes satisfy a functional central limit theorem 
(FCLT), that is, 

(2.1) Am. {D, Ml) as n oo, 
where =^ denotes convergence in distribution, 

(2.2) An{t) = C-\An{t)-Xnt), t>0, 

and {c„ : n > 1} is a sequence of positive numbers such that c,„ — )■ oo, n/c„ — 
oo and y/n/cn — )■ as n — )• oo. The canonical example is c„ = for 1 < 
a < 2. For background, see [1, 2] and [8]. 

As a consequence of this scaling, the arrival process satisfies the customary 
functional weak law of large numbers (FWLLN), that is. An =^ Xe in D as 
n — )■ oo where An{t) = n~^An{t) and e{t) = t for each t > 0. 
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When A „ is a renewal process for each n, the hmit process A = {j4(t) : t > 
0} win be a Levy process (have stationary and independent increments) . The 
Hmit then is naturally related to the FCLT for the sums of interarrival times, 
using the continuous mapping theorem and related arguments, associated 
with the inverse map together with centering (see Sections 7.3 and 13.7 of 
[8], especially Theorem 7.3.2 and Corollaries 7.3.2 and 7.3.3). 

The usual definition of the quality-and-efficiency-driven (QED) regime, 
leading to the celebrated square-root staffing rule, needs to be modified. For 
a modified QED regime, we assume that 

(2.3) c-^n{l- pn)^ /3, -oo</3<oo, 

where /?„ = X/nfi is the traffic intensity. With the heavier scaling here, the 
safety factor has to be greater: 1 — /9„ ~ f^Cn/n as n — y oo where ~ bn 
means that Un/bn — )■ 1 as n — )• oo. That implies a larger safety factor, because 
Cn/n goes to more slowly than 1/^/n. 

Let Qn = {Qn{t) :i > 0} be the queue-length process where Qn{t) is the 
number of customers in model n at time t. Define the scaled queue- length 
processes Qn = {Qn{t) : i > 0} and Qn = {Qn{t) : t > 0} by 

(2.4) Qn{t) = n-^Qn{t), Qn{t) = c-\Qn{t)-n), t>0. 

We can apply Theorem 1.1 above to establish a FCLT for Qn- The follow- 
ing theorem is an analog of Theorem 7.1 of [4] for the M/M/n + M model, 
and Section 7.3 of [4], which extends it to the G/M/n + M model, all with 
conventional QED many-server heavy-traffic scaling. We use Theorem 1.1 
with h[w) = — fi{w AO) — 9{w\/ 0) for all tt; G M where wiAw2 = m.m{wi,W2} ■ 



Theorem 2.1 (FCLT in the modified QED regime). Consider the model 
G/M/n + M in the modified QED regime (2.1)-(2.4)- If there is a random 
variable Q{0) such that QniO) =^ Q(0) as ?i— )• oo, then 

Qn^u} and Qn^Qin{D,Mi) as n—T-oo, 

where uj(t) = 1, t > 0, and Q = {Q{t) :t>0} is defined by the integral repre- 
sentation 

Q{t) = Q(0) - + A{t) - [ (/i(Q(s) A 0) + e{Q{s) V 0)) ds, t > 0. 

Jo 

Proof. As reviewed in [4], we first characterize the process Qn via the 
integral equation 

Qn{t) = Qn{0)+An{t)-s(fi [ (Q„(s)An)ds 

(2.5) 



-l{b j^{Qn{s)-n)+ d^, t>0. 



CONTINUITY IN THE Mi TOPOLOGY 5 

where (w)'^ = max{u;,0} and the processes S and L are independent rate-1 
Poisson processes (see Lemma 2.1 of [4]). By the definition of Qn in (2-4) 
and the integral equation in (2.5), we have 

Quit) = Qn(0) + i„(t) - Sn{t) - Ln{t) - ^lC-^n{l - Pn)t 

(2.6) 

{|I{Qn{s)^Q)+eQn{s) + )ds, 



where the processes Sn and !/.„ are defined by 

5„ = c^^ ^S'n ^/j, j {Qn{s) /\n)ds^ '^'^)'^^ 



As in Section 7.1 of [4], the processes Sn and L.„ are square integrable 
martingales with respect to the filtration F.„ = {J^nit) : t > 0} where 



Tn{t)=(jiySi^fi J {Qn{u) An)duj,Li^fi J {Qn{u)-nrduj:0<s<t 

Va{Qrr{0),An{s):s>0)VAf 

with M being the collection of all null sets. We cope with the general non- 
Markovian arrival process by putting the entire arrival process in the fil- 
tration. The predictable quadratic variation processes (Sn) and (L^) are 
defined by 

,4 

(5„)(i) = ^ / {Qn{s)/\l)ds, 
Cn Jo 
,4 

{Ln){t) = '^ {Qn{s)-l)^ds, t>0. 
Cn Jo 

By Lemmas 3.3, 5.8 and 6.2 of [4], the sequence of processes {Qn :n>l} is 
stochastically bounded in D. Applying the FWLLN for stochastic bounded 
sequences of processes in D in Lemma 5.9 of [4], we obtain the FWLLN: 
Qn =^ CO in D where ui{t) = 1, t > 0. Then by the continuous mapping theo- 
rem applied to the function (p: D ^ defined by 

ct){x){t)=(n [ {x{s)M)ds,e [ (x(s) - l)+(fsY t>0, 



and the assumptions on the scaling constants c„, we obtain {{Sn), {Ln)) =^ 
(?/, Vj) in D'^ as n — )• 00 where i]{t) = for all t > 0. By the martingale FCLT 
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(Theorem 7.1 of [2] and Section 8 of [4]), we obtain {Sn,Ln) =^ {fliV)- So we 
have the joint convergence 



(Qn(0), i„„ 5„, in) (Q(0), i, V, V) in M X D^, 



where is endowed with the product topology associated with the Mi 
topology on D. 

As noted before, the integral representation for Qn in (2.6) corresponds to 
(1.1) with function h{w) = —fi{w AO) — 9{w\/0) for all w S M. By continuous 
mapping theorem with the addition operation and the mapping in Theorem 
1.1, together with the convergence of the processes ((5.„(0), 5„, L„), we 
obtain the desired limit Qn =^ Q in {D,Mi). □ 

3. Proof of Theorem 1.1. We use the following elementary lemma, which 
we state without proof. 

Lemma 3.1 (Jump-coincidence). Given that y is the image of the map 
ip{x) defined in (1-1), the locations and sizes of the jumps of x and y must 
coincide. 

Proof of Theorem 1.1. Given that dM-i{xn,x) -)• 0, we let (itn,r„) 
and (u,r) be parametric representations of x„ and x, constructed as in The- 
orem 1.2, so that r and r„ are absolutely continuous with respect to Lebesgue 
measure with the properties in (1.4). Given these properties, we can follow 
the proof of Theorem 4.1 of [4] for the Ji topology, making appropriate 
modifications to cope with the Mi topology. 

The task is to construct associated parametric representations {uy^ , Vy^ ) 
and {uy,ry) for y„ and y. By the jump-coincidence lemma. Lemma 3.1, we 
can let ry = r and r^^ = r„ for all n. Then the desired convergence for the 
time components of the parametric representations follows from the assumed 
convergence Xn^ x: \\ry^ — ^y\\ = W^n — ^ as n ^ oo. Having specified 
the time components of the parametric representations of y and we must 
have 



where Disc(y) = {t: \y{t) — y{t—)\ > 0} and Disc(?/)'^ is the complement. 
Now suppose that s is such that r{s) S Disc(y)'^. Then from (1.1), 



Uy{s) = y{r{s)) 

UyAs) = Vnirnis)) 



if r(s) G Disc(y)'^, 

if r.„(s) G Disc(y„)'',n > 1 




■r(s) 



(3.1) 




u{s) + / h{y{r{w)))r' (w) dw, 
Jo 
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where the second hne follows by making the change of variables r(w) = z, 
so that r'{w)dw = dz (e.g., see Problem 13 on page 107 of [6]). 

In fact, because of Lemma 3.1, we can extend the representation in (3.1) 
to all s by simply letting 

Uy{s)=u{s)+ / h{y{r{w)))r' (w) dw, < s < 1. 
Jo 

Now observe that y{r{s))r' (s) = Uy{s)r'{s) for almost all s with respect to 
Lebesgue measure because r'(s) = whenever y{r{s)) ^ Uy(s). Hence we can 
write 



Uy{s) = u{s)+ / h{uy{w))r' (w) dw, 



0<s< 1. 



Applying the same reasoning to Uy^ , we can write 

UyAs)=Un{s) + [ h{Uy„{w))4,{w) dw , 



0< s< 1. 







Using these representations and the notation r]n = — we can then 
write 



\Uy„{s) - Uy{s)\ 

< \Unis) - u{s)\ + 



h{uy^{w))r'^{w)dw — / h{uy{w))r' (w) dw 



<rin + 



KuyAw))r'n{w) dw - / h{uy{w))r'^{w) dw 



(3.2) 



h{uy{w))r'^{w) dw — / h{uy{w))r' {w) dw 



<Vn + \\r'n\\c I \uy„{w) - Uy{w)\dw ■ 



\r'n{w) — r'{w) \ dw 



< {Vn+\\Hy)\\\K-r'\\i) + \\r'Jc j \uy^{w) - Uy{w)\ dw 
<(^n+||%)||||r;-r'||i)ell<ll-, 



where ||/i(y)|| = supQ<^<'p|/i(y(t))|, and we use the fact that h is Lipschitz 
with Lipschitz constant c together with (1.4) in the third line and apply 
Gronwall's inequality, as in Lemma 4.1 of [4], in the final line. Combining 
(1.4) and (3.2), we obtain 



\ny„-Uy\\<{7in + \\h{y)\\\K-r'h)e\\^'«\\^^0 



as n — )• oo. 



□ 
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4. Proof of Theorem 1.2. We break up the proof into parts presented in 
the following subsections. First, in Section 4.1 we establish some bounds on 
the maximum jump function and the uniform (or supremum) norm 

||x|| = sup 

0<t<T 

In Section 4.2 we show that the time component r of the parametric repre- 
sentation (n,r) of the limit function x can have the desired representation 
in Theorem 1.2. We then define the associated function u in the parametric 
representation (n, r) of x as required, using linear interpolation where there 
is freedom. 

We construct the desired parametric representations of Xn in Sections 
4.3-4.7. In Section 4.3 we partition the domain [0, 1] into finitely many 
subintervals, of which there are three kinds. We then carry out the proof 
for each of the three kinds. We obtain this finite number by considering the 
finite number of discontinuities of x exceeding some small £i . The first kind 
of subinterval corresponds to the flat spots in r associated with the large 
(of size bigger than ei) discontinuities in x. The second kind of subinterval 
corresponds to short connecting open subintervals between the closed subin- 
tervals with large jumps and the closed subintervals with no large jumps. 
The third kind of subinterval corresponds to subintervals where there are 
no large jumps, but there may be (even infinitely many) small jumps. We 
construct the new parametric representations for these three kinds of inter- 
vals in Sections 4.5-4.7. In each case we show the convergence as required 
for the metric in (1.3) with the extra properties in (1.4). In Section 4.4 we 
show how to construct the spatial part of the parametric representations of 

Xn ■ 

4.1. Bounds on the uniform norm. For treating the closed intervals with- 
out jumps exceeding a small threshold, we apply some bounds on the maximum- 
jump function and the uniform norm, which may be of independent interest, 
so we establish them first. We will use the uniform norm for real- valued func- 
tions with different domains; the desired domain should be clear from the 
context. Normally, when we write ||3;|| for x € D, the domain is understood 
to be [0,T], but when we write \\u\\ and ||r|| for a parametric representation 
of X, the domain is usually understood to be [0, 1]. However, we will also con- 
sider the uniform and Li norms over subintervals; the relevant subinterval 
should be clear from the context. 

A key role is played by the maximum-(absolute)-jump function. Let 

(4.1) Jmax(x) = snp{\x{t) - x{t-)\ :0<t<T}. 

(The supremum in (4.1) is always attained because, for any e > 0, \x{t) — 
x{t—)\ > £ for only finitely many t in [0,T].) 
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Lemma 4.1. Ifxn^xin {D,Mi), then 

limSUp{Jmax(a;n)} < Jma.^{x). 

Note that we need not have equahty in Lemma 4.1, because the functions 
Xn could have smaher jumps. Indeed, the functions Xn might be continuous, 
in which case the hmsup is zero. 

Proof of Lemma 4.1. We will show that a subsequence of the locations 
and sizes of the maximum jumps of x„ necessarily converge to a limit for x 
which provides a lower bound for the maximum jump of x. We will exploit 

compactness to obtain convergent subsequences. Let {on^} = {a^^ :n > 1} 

denote a subsequence of the sequence {an} = {a„:n > 1}, and let {an } 

denote a subsequence of the subsequence {an^}, and so forth. 

Given dMx{xn-,x) — )■ 0, we can choose (u,r) G n(x) and {un,rn) G n(x„) 
for n > 1, such that ||n„ — u\\ V ||r„ — r|| — )• by Theorem 12.5.1(i) in [8]. For 
each n > 1, let si^n and S2,n be points in [0, 1] such that 

dn = \Unis2,n) " Unisi,n)\ = JmaA^n), n>l. 

[Recall that the supremum in (4.1) is attained.] Choose a subsequence {dn^} 
of {dn} such that 

d^^^ — )• lim sup { Jniax (xn)} asn— )-oo, 

n— )-oo 

which is always possible by the definition of the lim sup. 

Let {s^^^} and {s^l} be the associated subsequences of the original se- 
quences {si,n} and {s2,n} yielding the sequence {dn^}. Let tn^ = rn\si^n) = 
i^n\s2,n) be the associated flat spots. Since tn^ is an element of the compact 

(2) 

set [0,T], there exists t G [0, T] and a subsequence {tn } of the subsequence 
{tn^} such that tn^ — )• t as n — >■ oo. 

(2) (2) 

Let {s\ and be the associated subsequences of the subsequences 

{^i^ni ^^"^ {•^2^1} corresponding to {tn^}. We can thus find points si,S2 G 

(S) (3) 

[0,1] and further subsequences {^i „} and {s2n} of these subsequences so 

(3) (3) 

that s\ — )• si and n ~^ •^2 as n — )■ oo. However, by the uniform convergence 
{un,rn) to {u,r), we have the associated limits 

ri?Hsfl)^r{si)=t and rl^^s^^^) ^ r{s2) = t asn^oo, 
ul^Hs^iJi) ^ u{si) and u^f'' (s2^^) — ^ ii(s2) asn— ^oo, 

d^n^ = |n(fH4'l) - ^^^nHsfh ^ limsup J^ax(Xn). 
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Since r{si) = r{s2) = t, we necessarily have 

JmaAx) > \xit) -x{t-)\ > \U{S2) - u{si)\ = lim SUp Jmax(a^n.) • □ 

n—^oo 

Next we introduce several oscillation functions. As in (2.5) (on page 393, 

(3.1) on page 394, (4.4) on page 402 and (5.1) on page 404 of [8]), let 

u{x,A)= sup {\x{ui) — x{u2)\}, 
u{x,6)= sup {u{x,[t,t + S))}, 

0<t<T~S 

(4.2) \c-[a,b]\= sup {\c- {pa + il-p)b)\}, 

0<p<l 

Ws{x,t,5)= sup {\x{t2) -[x{ti),x{t3)]\}, 

0V{t-5)<ti<t2<t3<{t+5)AT 

■Ws{x,6)= sup Ws{x,t,6). 

0<t<T 

In our proof of Theorem 1.2, we exploit part (b) of Lemma 4.2 below. 

Lemma 4.2 (Maximum-jump bound on the uniform norm), (a) Suppose 
x,Xn^ D with {u,r) G n(x) and {un,rn) S n(a;„). Then 

(4.3) \\Xn- x\\<Ws{Xn,\\rn '"ID 2</max(^n) ~l~ 'Anax(^) ~l~ W'^n 

(b) Ifxn^x in {D,Mi) then 

(4.4) limsup{||x„ - x\\} < 3Jmax(a;)- 

n— >-oo 

(c) //, in addition, continuous, then 

limSUp{||Xn - xll} < Jmax(a;)- 

n— >oo 

Proof. By the triangle inequality, 

\\xn — x\\ = \\Xn o r — X o r\\ 

< \\Xn O r — Xn O + \\Xn o Tn — Un\\ + \\Un — u\\ + ||n — X O r||. 

We bound the first term on the second line by observing that 

\\Xn or-XnOrn\\< l^{Xn, - r\\) < Ws{Xn, \\rn - r\\) + Jmax(a;n,), 

using the definitions in (4.1) and (4.2). Since o r„ — Un\\ = ^max(a^n) and 
||n — xo r|| = Jmax(2;), that explains (4.3). 
We now turn to part (b). First, we note that 

(4.5) limlimsupt(;s(x„,e) = 0, 

£4-0 n— >oo 
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by Theorem 12.5.1(iv) of [8]. Since — )• a; as n — ?• oo, we can select para- 
metric representations such that — "uHvUrn — r||— )'0 as n— )-oo. To- 
gether with the hmit in (4.5), that imphes that Ws{xn, \\rn ~ A\) ^ 0- Since 
II — u\\ -^0 as n— )-oo, the hmit in (4.4) follows from Lemma 4.1. Finally, 
part (c) follows easily from parts (a) and (b) because Jmax(2;n) = when 
x„ G C. □ 

4.2. Constructing the parametric representation of the limit function. We 
start by constructing a special parametric representation (u, r) of the limit 
function x. 

Lemma 4.3 (Choice of r in the parametric representation oi x). For any 
X G L'([0,T],M), there exists a parametric representation {u,r) of x such that 
r is absolutely continuous with respect to Lebesgue measure, having deriva- 
tive r' almost everywhere, satisfying \\r'\\ < 2T. Moreover, this function r 
can serve to build the parametric representation {u,r) of x needed to es- 
tablish convergence dMi{xn,x) — )■ for any sequence {xn'-n > 1} for which 
convergence holds. 

Proof. By Theorem 12.5.1(i) of [8], there is total freedom in the choice 
of the parametric representation {u,r) of the limit function x. We can start 
with any proper parametric representation of x, and if convergence Xn^ x 
holds, then it will be possible to find suitable parametric representations 
{un,rn) of Xn- So the construction we carry out for r and u are necessarily 
without loss of generality as far as establishing the convergence is concerned. 
However, we need to show that it is possible to find a parametric represen- 
tation of X with the additional structure. 

Given that r: [0,1] — >• [0,T] is onto, the maximum value of its derivative 
(if it exists) must be at least T. Indeed, when x is continuous, we can just let 
r'{s) = r, so that r{s) = Ts, < s < 1. However, there is no need for the Mi 
topology unless the limit x has at least one jump. So henceforth we assume 
that is the case. Then the parametric representation must have a fiat spot for 
each jump; that is, if x has a jump at t, by which we mean \x{t) — x{t—)\ > 
0, then there must exist an interval [si,S2] C [0,1] such that r(s) = t for 
s E [si,S2]; ''(s) < i for s < si, and r{s) > t for s > S2. To concisely express 
that, we write r~^{t) = [si, S2] where r~^{t) = {s G [0, 1] : r{s) = t}. We now 
show that r can be constructed by combining linear pieces such that each 
piece either has slope (is a fiat spot) or has slope 2T. We let r have a flat 
spot at t for each t G Disc(3;). It is possible that the interval endpoint T is a 
discontinuity point of x. Whether or not T G Disc(x), we include a flat spot 
at T in r. Otherwise, r has no flat spots at t £ Disc(3;)'^. We guarantee that 

< ^('^2) - ^-(51) ^ for < si < S2 < 1 and n > 1. 

S2 - si 
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The construction is elementary if x has only finitely many discontinuities, 
but the number of discontinuities can be countably infinite, even dense in 
[0,1]. Thus to carry out the construction, we initially order all the discon- 
tinuity points of X in order of decreasing size of the jumps so that ti is the 
location of the largest jump while t2 is the location of the second largest 
jump, and so forth. We can break ties arbitrarily. To be definite, suppose 
that ties are broken by taking the discontinuities in order of their time value; 
for example, if the jumps at the points ti and t2 are the same size, then we 
order them so that ti < ^2 • We assign a flat portion to r of length fj for the 
jth discontinuity. We choose these lengths fj such that their sum is 1/2. 
That leaves a total length of 1 /2 in the domain [0, 1] to be the support of 
the increase of r. Wherever r can have an interval without discontinuities, 
we let r increase at slope 2T. Thus ||r'|| = 2T. 

We start by letting the successive lengths of the flat spots, when first 
introduced in the construction process, satisfy the inequalities 

oo 

(4.6) > ^ ^ 1; 

i=j+i 

that is, the length of each flat spot exceeds the sum of the lengths of all 
remaining flat spots. (For example, that can be achieved by letting fj = 
2~^-^, j > 2, and fi = 5/12.) Inequality (4.6) allows us to remove portions 
of an initially assigned flat spot, taking away the length of a new flat spot. 
Requirement (4.6) ensures that there is a positive length for each flat spot 
after all these subsequent changes. There will also be a flat spot at T, making 
the sum of all flat spots be 1/2. 

We construct our function r iteratively. In particular, we will construct a 
sequence of functions {r„ : n > 1} (not to be confused with the parametric 
representations of Xn discussed later) such that r as n — t- oo. Let A be 
Lebesgue measure on the interval [0,1]. For each n, r^(s) will equal either 
2T or for almost all s with respect to Lebesgue measure. We will further 
have 

oo 

(4.7) A({s:r;(s)/r'(s)}) = 2 ^ /,^0 asn^oo. 

i=n+l 

As a consequence, r' will also equal either 2T or for almost all s. An 
example of the first four steps of the construction to be described, yielding 
r4, is shown in Figure 1. 

Specifically, to construct r, we start with ro(s) = 2Ts, < s < 1/2, and 
ro(s) = T, 1/2 < s < 1. Flat spots with value T will be present in r„ for all n 
and in r. To construct the next function ri, we append a flat piece of length 
/i to the initial function rg, extending out to the right of length fi at ti on 
the y axis. The flat portion starts at si = ti/2T and extends to si + /i. We 
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time component of the parametric representation 




s= 1/2 



Fig. 1. Constructing the time component r : [0, 1] [0,T] of the parametric representa- 
tion (u, r) of the graph Fx of a function x. We show the function r^ which is the result of 
the first four steps of the construction. The functions r„ are all piecewise linear, with r„ 
containing 2n + 1 pieces, each either of slope or slope 2T. 

then add an increasing portion with slope 2T, starting at the right endpoint 
si + /i. The new function is the lower envelope of these pieces. 

We continue in this way to construct a sequence {r„ :n > 1} of functions 
which decrease toward the limit r. The construction of r„ given r„_i pro- 
ceeds in the same way when tn > tj for all j , 1 < j < n — 1 ; that is the case 
for t2 in Figure 1. Otherwise, like in Figure 1, the linear piece added at 
Sn + fn will hit r„_i on the y axis at tn = min{tj -.tj > tn,l < j < n — 1}. 
For in Figure 1, that occurs for ^3 = ti. Then r„ coincides with r„_i 
except for two new linear pieces: first, the flat spot mapping [smSn + fn] 
into tn and, second, the linear piece with slope 2T, r„(s„ + fn) = tn and 

'^ni.^n ~l~ fn ~l~ (^n ^n)/2T') — tn- 

Note that the functions r„ are all piecewise linear, with r„ containing 
2n + 1 pieces, each either of slope or slope 2T. Thus is absolutely 
continuous with respect to Lebesgue measure for each n. For all n and almost 
all s, r'^{s) G {0,2T}. In addition, r„(s) > r„+i(s) > 0. By this construction, 
the flat portions at any level tj may decrease as flat portions are added 
for t < tj, as is the case for ti in Figure 1 when we add the flat portion 
of length /s corresponding to t^. If we add a flat spot of length fj at step 
j, then the total length of the change in the derivative is 2fj, which gives 
(4.7). By condition (4.6), there always will be a remaining flat portion at 
tj whenever one was initially inserted there. Since, r„ is strictly decreasing. 
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there necessarily exists a function r such that r„(s) — )■ r as n — t- oo for all s. 
In addition, we have convergence of the derivatives, as in (4.7). 

In general, the values r'^{s) could oscillate among the two values 2T and 
infinitely often, but 

Q^ rn{s2)-rn{si) for 0<si<S2<l andn>l. 

S2 - Si 

Indeed, this construction guarantees (4.7). As a consequence, the sequence 
{r'n : > 1} converges to r' in Li, that is, as n — )• oo, 

Wr'n - r'Wi = K{s) - r'{s)\ ds < 2T\{{s : r'^{s) / r'{s)}) ^ 0. 
Jo 

We now specify the space component u of the parametric representation 
of X. Recah that Disc(x) = {t £ [0,T] : \x(t) - x{t-)\ > 0}. We must have 
u{s) = x{r{s)) for each s with r(s) = t £ Disc(x)'^ and r(l) = T. [We could 
have T G Disc(x), but we necessarily have G Disc(a;)'^.] It remains to specify 
ti at s if r{s) = t G Disc(x) (with t^T). Whenever [si, S2] = r~^{t) for S2 > 
si, let u{si) = u{si — ) and u{s2) = n(s2+), where the left and right limits 
are over s such that r(s) = t G Disc(x)'^. Since x has left and right limits 
everywhere, these limits are well defined, with u{si) = x{t—) and ^(^2) = 
x{t). Then let the remainder of u be defined by linear interpolation, that is, 

u{s) = pu{si) + (1 - p)u{s2) = px{t-) + (1 - p)x{t) 

for all s = ps2 + {I — p)si, <p <1. Thus we have constructed the desired 
parametric representation {u,r). □ 

Given the constructed parametric representation {u,r) of x, and given 
that x„ — )• X in (D,Mi), let (ii„,r„) be a parametric representation of Xn, 
n > 1, such that \\un — u\\ V ||r„ — r|| — )■ as n — )■ 00, which must exist, 
by Theorem 12.5.1 of [8]. Our goal, then, is to construct new parametric 
representations {un,fn) of x„ such that f„ has the properties in Theorem 
1.2, including (1.4), and ||tt„ — u\\ V ||f„ — r|| — as re — )• 00. 

Our goal can be expressed as showing that there exists re** = re**(e, u, r, {xn'-, 
n > 1}) for any specified e > 0, {u,r) G n(x) as constructed above and 
{xn;n > 1} where dMi{xn,x) — )■ 0, such that, for all re > re**, there exist 
parametric representations {un,^n) of Xn, such that f„ has the properties in 
Theorem 1.2, including (1.4), and 

(4.8) ||m„ — u|| V ||f„ — r|| < e. 

Indeed, we will construct such an re**, with the final specification being 
(4.30). There will be several steps. 
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4.3. Constructing the finite partition of the domain. Here is where we 
start: We are given the hmit x and the sequence {a;„ :n > 1} with x„ — t- x 
as n — 7- oo in Z)([0,r],M) endowed with the Mi topology. Let (u,r) be the 
parametric representation of x constructed in Section 4.2. Let {un,rn) be 
the parametric representations of x„ such that ||u„ — m|| V ||r„ — r|| — )• as 
n — )• oo. Let e > be given. 

Choose El <e/9. [The reason for this inequahty is explained at the very 
end, in (4.31)]. To simplify the presentation, assume that x has no dis- 
continuity at T. It is not difficult to treat that case too, but it is slightly 
different. Let Disc(x,ei) = {t € [0,r] : \x{t) — x{t—)\ > ei}. Let m = m{ei) = 
|Disc(x,ei)|, the cardinality of the set Disc(x,ei), which is necessarily finite 
(see Theorem 12.2.1 of [8]). Let these jump times be labelled, so that 

= to < ti < ■ ■ ■ < tm < tm+i = T. 

Note that this labelling is different from the labelling used in Section 4.2. 

We next introduce £2 and n{e2) so that — n|| V ||r„ — r|| < £2 for all 
n > n{e2). With that in mind, we now choose 62 such that the following four 
properties are satisfied: 

(i) 0<e2 <ei, 

(4.9) tj + 2e2<tj+i-2£2, l<i<m-l, 

< ti - 2^2, and tm + 2^2 < T] 
(ii) 

(4.10) i/(x, +362)) < ei/2 and v{x,[tj - ■ie2,tj)) < ei/2 

for 1 < j < m, where i' is the modulus of continuity defined in (4.2); 

(iii) There exist time points S2j_i{£2) and S2j(e2), ^ ^ j ^rn, such that 

(4.11) r~\tj-2e2) = {s^._,{e2)} and r~\t, + 282) = {s+ {62)} 

for 1 < j < m, that is, s^__]^(e2) and S2j{s2) are points of increase (not in 
flat spots) for the function r (tj it 2^2 are continuity points of x); 
(iv) 

(4.12) S2j-i-s2j_i{e2)<-^ and s^-{e2) - S2j < 

for 1 < j < m, So = and S2m+i = 1, where r~^{tj) = [s2j-i, S2j], 1 < i < m. 

The second property holds for all sufficiently small £2 because of right 
continuity and the existence of left limits everywhere. As a consequence of 
property (4.9), 

0< 5^(62), 4i(e2) <s^j+i(e2), l<j<m-l, s+„(e2)<l. 
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time component of the parametric representation 




Fig. 2. An example of the partition of the domain of the time component r: [0, 1] — [0,r] 
of the parametric representation {u, r) of the graph Vx of a function x into 4m + 1 dis- 
joint subintervals, associated with the m large jumps. In this example, there are m = 2 
large discontinuities in [0,r], so that the domain is partitioned into 4m + 1 = 9 disjoint 
subintervals. 

We can obtain all these properties because we are imposing only finitely 
many requirements. For condition (4.11), we use the fact that x has at most 
countably many discontinuities. 

As a consequence of the construction above, we have partitioned the do- 
main [0, 1] of {u,r) into 4m + 1 disjoint subintervals: 

[s2j~l,S2j], l<j<m, 

(%_i(e2),S2i-i) and (s2i,S2j(^2)), 1<J<?TI, 

[0,sf(e2)],[4m(e2),l] and [s^j(e2), s^j+i(e2)], l<j<m-l. 

The case of m = 2 is shown in Figure 2. The m closed subintervals [s2j^i,S2j] , 
1 < j < m correspond to the m large (of size bigger than ei) discontinu- 
ities of X, and the associated m flat spots of r. Corresponding to each of 
these m subintervals, we have two connecting open subintervals. Thus, over- 
all, there are 2m connecting subintervals of the form (s^_^(e2)) ■52j-i) and 
{s2j,S2j{£:2)), 1 < i < Finally, we have m + 1 closed subintervals in which 
there are no jumps exceeding ei: [0,sj"(e2)], [s2m(^2),l] and [s2j(^2), S2j+i('^2)], 

1 < i < ""1- — 1 • 

Before proceeding, we explain the strategy of our proof and the de- 
tails in properties (4.9)-(4.12): the use of 2^2 in (4.9), 3^2 in (4.10) and 
£i/6mT in (4.12) (see Section 4.6 for more details). We will let be a 
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minor modification of r„ and r, respectively, in tlie m large-jump closed 
subintervals [s2j-i,S2j], and in tlie m + 1 small-jump closed subintervals 
[s2j{^2) , S2j_^_i{£2)]- We use the remaining 2m connecting open subintervals 
to make adjustments ensuring that f„ is nondecreasing with the desired 
properties. Condition (4.10) ensures that these small open subintervals are 
manageable. 

Toward that end, we let 

rnis2j-l) = rn{s2j-l) > r{s2j-l) -£2= tj - £2, 

(4.13) 

rn{s2j) = rn{s2j) < r{s2j) +82= tj + £2 

for all j and n>n{£2), but 

rn{s2j-i{e2)) = r(s^-_i(e2)) = tj - 2e2, 

(4.14) 

r„(4j(e2)) = r{stj{£2)) = tj + 2£2 

for all j and n > n{£2). Note that f„ is defined in terms of r„ in the first 
two cases and in terms of r for the second two. 

We have used 2e2 in (4.9) to ensure that fn{s2j-i) > ^n('52j-i(^2)) and 
fn{s2j) < rn{s2j{£2)) for n > n{£2). We will typically need strict inequality 
in order to construct f„ properly in the subintervals {s2j^i{£2), S2j-i) and 
is2j , s^j{£2)) for 1 < j < m. 

We use 3^2 hi (4-10) to control the oscillations of Un in the small open 
intervals (s^_i(e2); S2j-i) and (s2j, S2j(e2)) for l<j<m. We intend to let 
Un{s) = Un{fn{s)). Hencc, we apply (4.10) and (4.14) to obtain, first, 

rn(s2i-l(^2)) = r{s2j_i{e2)) = tj - 282 > tj - 3£2, 

rn{s^j(.£2)) = r{s^j{£2)) = tj + 2£2 < tj + 3e2 

and, second, 

tj -£2> ?'n(s27-l(^2)) > tj - 3£2, 

(4.15) 

tj +£2< rnis2j_i{£2)) < tj + 3^2 

for all j and n > n{£2). 

Finally, we use £i/6mT in (4.12) to make these 2m open subintervals 
(■^27-1(^2)1 •S2j-i) and {s2j, S2j{£2)) for I <j <m so short that they all can 
only contribute £1 to the Li norm Wr'^ — r'\\i when < ||r'|| < 2T and < 
llr^ll < 3T [see (4.26) in Section 4.6]. 
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4.4. The spatial component of the parametric representations. In this 
subsection we show how to construct Un, the spatial component of the new 
parametric representation of Xn- As a basis for having a simple construction, 
we require that f .„, the temporal part of the spatial representation of Xn, have 
its flat spots be in one-to-one correspondence with the flat spots of r„. The 
given function r„ must have a flat spot corresponding to each discontinuity 
point of Xn, but it could have additional flat spots (which could be removed, 
but we do not do that). We take r„ to be as given, and set as a requirement 
for f„ that its flat spots correspond to those of r„. Specifically, We define 
f„ so that 

fn =rnO(t>n, 

where 4>n is an increasing homeomorphism of the domain [0, 1] . We will be 
constructing f„ and in the remaining sections. To achieve the desired 
order, as specified for (pn, we will do the construction in each case so that 
the flat spots in f„ correspond to those of r„ and appear in the same order 
(just as was done for r in Section 4.2). Given that direction, here we are 
showing how to construct Un given f„ and 

As part of our construction, we directly relate the flat spots: If r~^(t) = 
[si, S2] where si < S2, then f~^{t) = [si, §2] where si < S2, and (j)n{si) = Sj for 
i = l,2. Moreover, for each flat spot, we define (pn within such a subinterval 
by linear interpolation; that is, we let 

(4.16) 0n(ps2 + (1 -p)si) =ps2 + il-p)si- 

Finally, we construct and (pn so that the order of the flat spots is pre- 
served. (To do that, we can order the flat spots of the given rn in order of 
their length. We can do the construction, as in Section 4.2, preserving the 
desired order at each step.) Whenever, r^^(f„(s)) is a one-point set, we have 

(4.17) Ms) = r-HUs))- 

Given the homeomorphism (pn described above, it is elementary to define 
the function Un'- We let 

(4.18) Unis) = {UnO (j)n){s) = Un{(l)n{s)), < .S < 1. 

Whenever r~^(f„(s)) is a one-point set, we have 

Un{s) = {Un O (t)n){s) = {Un o O f„)(s) 

= {Xn o r„ o r~^ o f„)(s) = {Xn o f„)(s) = Xnifn{s)). 

With Un defined in this way, we will want to bound — We will do 
that by applying the triangle inequality: 

\\Un — u\\ = \\Un ° 4>n — ""11 < \\Un ° (pn — U ° (pnW + \\u ° 4>n — u\\ 

(4.19) 

< \\Un — u\\ + \\u o <pn — u\\ 

and work to control the second term on the right, for example, using (4.10). 
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4.5. Case 1: The flat spots corresponding to the large jumps. In this 
subsection we consider the m closed subintervals [s2j-i-,S2j\, 1 < J < m, 
corresponding to the m identified "large" jumps of x. We have r[s) = tj for 
all s G [s2j-i,S2j]- Since r is constant on these intervals, these intervals are 
easiest to treat. We are given r„ such that ||?^n — '^ll — '•0 as n— t-oo. However, 
these functions need not satisfy the smoothness properties of Theorem 

1.2. Thus we will construct new parametric representations (uni^n) over 
this subinterval. We will let f„(s) = rn(s) at the end points S2j-i and S2j 
of the interval [s2j-i,S2j] and also at finitely many points sj^k within the 
subinterval [s2j~i, S2j], but we will make a new definition at other points s, 
using a minor modification of the construction in Section 4.2. The end points 
S2j-i and S2j have been specified in the process of choosing 82 in Section 

4.3. We choose the internal points Sj^^, < k < nj, so that ^ < Sj,fc+i for 
all k and the constructed subintervals [sj^k,Sj^k+i] are of equal width 

(4.20) Wj = Sj^k+i-Sj^k = — — , 

no- 
where Uj will be specified below in (4.22). We will let fn{sj^k) = 'i'n{sj,k) and 
Un{sj,k) =Unisj,k) for ah j and k. 

Since we are redefining our parametric representations within these small 
subintervals, we want to construct the widths wj sufficiently small so that 
the fluctuation of u within the interval [sj^k, Sj^k+i] is suitably small. Recall 
that u has been defined by linear interpolation. Thus 

\u{Sj,k+l) - u{Sj,k)\ = . 

Uj 

Let Uj be chosen for each j so that 

(4.21) i/(n; [sj^k, sj^k+i]) = \u{sj^k+i) - u{sj^k)\ < £2/2; 
that is, let 

(4.22) nj = \2\x{tj)-x{tj-)\/e2], 

where [t] is the least integer greater than or equal to t. Then the subinterval 
width Wj is obtained by combining (4.20) and (4.22). 

By focusing on u, which depends on x, we have specified nj and Wj. It 
remains to ensure that the slope of f.„ need not exceed 2T. We let 

(4.23) £3 = min{e2 /2m, T mm{wj :1 < j < m}/2}. 

[We explain this choice after (4.24) and (4.25) below.] Given the above, we 
choose n* = n*{es), so that dMi{xn,x) < £3 for n > n*{es). (We will need to 
impose yet another constraint on n [see (4.30)].) 

We now can construct our desired parametric representations {un,^n) of 
Xn such that \\un — u\\ V ||f„ — r|| < 62 for n>n* = n*(e3), specified above. 
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where the uniform norm is restricted to the intervals [s2j-~i,S2j], 1 < i < 
As indicated above, we let fn{sj^k) = i'n{sj,k) for all j, k and n. Our idea is 
to use the construction in Section 4.2 within each subinterval [sj^fc) to 
construct f„ at other s, but we modify the construction slightly. In Section 
4.2, we constructed r to have a flat spot corresponding to each discontinuity 
point t of X. However, as indicated in Section 4.4, here we construct 
so that it has a flat spot for each flat spot in r„ (which may or may not 
correspond to the same jump point of Xn)- Since r„ has at most countably 
many flat spots in all, we can apply a variant of the construction in Section 
4.2 to them. We can order these flat spots of r„ according to their length in 
the domain. 

First, we let f„(s) =r„(sj,fc) for ah s £ [sj^k, Sj,k+i] if ''nlsj.fc) = ?^n(sj,fc+i)- 
Hence it suffices to consider subintervals for which rn{sj^k) < ^n(sj,fc+i)- We 
adjust the construction in Section 4.2 in the obvious way to allow for the 
different values at the endpoints in this context. In particular, here we have 

tj - £2 = r{s2j~l) - £2 < rn{s2j^l) < rn{Sj^k) 

< rniSj,k+l) < rn{s2j) < r{s2j) + £2 = tj + £2 

for all k and n > n(£2). Hence, by the definition of £3 in (4.23), for n > n*(£3), 
r-n(gj,fc+i) - rn(Sj,fc) ^ 2£3 ^ ^ 

Sj,k+1 — Sj^k ~ Wj ~ 

We initially let fnfl{{sj^k + •5j,fc+i)/2) = rn{sj^k+i), which makes the maxi- 
mum possible slope in the subinterval be 2T, just as in Section 4.2. We then 
carry out the same complete construction on each subinterval, getting f„^/ 
for Z > 1, and then the limit 'rn,oo = (on this subinterval [sj^k, Sj^k+i])- 

That construction gives us the function f„ which is absolutely continuous 
with respect to Lebesgue measure, having derivative almost everywhere, 
satisfying ||f^|| < 2T for n > n*(£3) where the uniform norm is restricted 
to the intervals [s2j-i) S2j], 1 < j < m. Moreover, since r is constant on the 
subinterval [s2j _ 1 , S2j ] , 



(4.24) 



2i fS2j 

\r'n{w) — r' (-01)1 dw = / r[^{w)dw 

2j-l J ^-23-1 



< 2£3(s2j - S2j-l) < 2£3 < — 

m 



for all n > n*(£3). The sum of these integrals over the m intervals [s2j-i, S2j] 
is then bounded above by £2. Thus — < £2 for all n > n*(£3) where 
the Li norm is restricted to the union of the intervals [s2j-i,S2j], 1 < J < m. 

Moreover, since dMi{xn,x) < £3 for n > n*(£3), ||r„ — r\\ < £3 for n > 
n*(£3), which implies — r|| < £3 for n > n*(£3) by the construction of 
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f„ above from r„, where the uniform norm is restricted to the intervals 
[s2j~i,S2j], l<i<m. 

We have akeady indicated how to construct Un in Section 4.4. Formulas 
(4.16)-(4.18) imply that n„ is defined by linear interpolation within each 
flat spot, just like u. Hence {un,rn) is a parametric representation of x.„ for 
each n, restricted to the given subinterval [s2j-i, S2j]- 

Finally, we need to bound — By (4.19) and (4.21), on the subin- 
terval [s2j-l,S2j], 

(4.25) \\un - u\\ < \\un - u\\ + max{ i/(n; [sj^k, Sj,k+i])} < £3 + tt < ^2 

k Z 

for n > n*(e3) by (4.21) and our assumption that n > n*(e3) [which implies 
that dMiixn,x) < £3 for those n and those subintervals] . Thus, we have 
established (4.8) for n > n*{es), where £3 is given in (4.23), which in turn 
depends on Wj (or, equivalently, nj) and £2, and where the uniform norm is 
restricted to the intervals [s2j-i, S2j], l<j<m. 

4.6. Case 2: The connecting open subintervals. We have introduced the 
connecting intervals (s2j_i(e2)) •S2j-i) and (s2j, s^(e2)) in order to provide a 
bridge between the closed intervals with large jumps and the closed intervals 
without large jumps. We have already let rn{s2j-i) = rn{s2j-i), fn{s2j) = 
rn{s2j), rn{s2j_i{e2)) = r{s2j_i{e2)) and f„(4j(e2)) = r (5^^.(62)) for all j 
and n > n{£2) in (4.13) and (4.14). These short open subintervals allows us 
to make the transition between these two different definitions. By condition 
(4.12), we have made these intervals short, so that their contribution to the 
total Li distance — can be controlled without carefully examining 
the derivatives. We only discuss a typical "lower" interval, (s^j_i(e2), ■52j-i), 
because the associated "upper" interval, {s2j, S2j{£2)), can be treated in 
essentially the same way. 

We already have observed that (4.9) ensures that f„(s2j-i) > f^n(s2j_i(^2))- 
We initially define f„ by linear interpolation between the established assign- 
ments at the interval endpoints. We then make asymptotically negligible 
adjustments by adding flat spots as necessary to account for the given flat 
spots of r„ (which includes all discontinuity points of Xn in this range). Note 
that f„ maps the subinterval [s2j_i{£2),S2j-i] into some interval 
We include flat spots in f„ with domain in [s2j_i{£2), S2j-i] to match all flat 
spots of r„ with values in [tnjjf^j]. 

Since ||r„ - r|| < £3 < £2/2 for n > n*(£3), < r(s2j-i) - r(s^^_-^(£2)) = 
2£2, f„(s2j-i) = r„(s2j-i) and f„(s2j_i(£2)) = r„(s2j„i(£2)) by the construc- 
tion, we must have f„(s2j_i) — ^n('S^„x(£2)) < (5/2)£2 and ||f„ — r\\ < 3£2 
for n > n*(£3) where the uniform norms are restricted to the connecting 



22 



G. PANG AND W. WHITT 



open subintervals. Moreover, before the addition of any flat pieces, over the 
open subinterval (s2j_i(^2)5 •52j-i), 

„ r„(s2j-i) -r„(s2,-_i(e2)) (5/2)62 
r = < — — 

S2j-1 - S2j_iie2) S2j-1 - S^j_i(e2) 

S2,^i-s^j_,{82) J 4 II'' II 2 

for n > n*(e3). Hence we can achieve ||f^|| < 3T over the subinterval in 
the first and the second cases for all n > n*(e3) even after adding small flat 
pieces corresponding to the flat spots of r„ and thus making the linear pieces 
slightly steeper. 

Since < ||r'|| < 2T and < ||f^|| < 3T for n > n*(e3), we have the crude 
bound ||f^ — r'll < 3T on the subinterval (s^_;^(e2)) ■52j-i) for n > n*{e3). 
Since the length of this subinterval is very short, with bound in (4.12), on 
this subinterval, we have 

(4.26) \rM - r'{w)\ dw < 3T{s2j-i - s^^^,{e2)) < ^ 

for n > n*(e3). Hence the total contribution to the overall Li norm ||f,^ — '"'Hi 
from all the 2m connecting open intervals is bounded above by ei . 

We have already indicated how to construct the spatial portions Un of 
the parametric representations of Xn in Section 4.4. However, we modify the 
bounding argument in (4.25). Since (4.14) and (4.15) hold for all n > n*(e3), 
and by construction r, r.„ and f„ are all nondecreasing over the interval 
(^27-1(^2), •S2j-i), we must have 

tj - 3^2 < r{<j)n{s)) < tj for 52^-1(^2) <s< S2j~i 
for n > n* (e^) . Hence we can apply the oscillation bound in (4.10) to obtain 

(4.27) \\uo(j)n - u\\ < v{x, [tj - 3e2,tj)) < ei/2 

for n>n*[£^), where the uniform norm is over the connecting open subin- 
tervals. Applying (4.25) and (4.27), we get 

II II ^1 

\\un - n|| < £3 + y < ei 

for > n*(e3), over the connecting open subintervals. 

4.7. Case 3: The closed intervals without large jumps. We now treat the 
m + 1 closed intervals [527(^2), •s^_,_]^ (£2)], < j < m, corresponding to the 
portions of x without any large jumps. Here we do not pay careful attention 
to the jumps of the limit function x. The idea is that the resulting error is 
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bounded above by a constant multiple of the size of the largest jump in the 
limit function x in this region, and is thus controlled by the fact that the 
size of each jump is at most ei, as stipulated in Section 4.3. 

The smoothness properties required by Theorem 1.2 are achieved by let- 
ting f„, the time portion of the parametric representation {un,fn) of x„, be 
a minor modification of r, which has already been shown to have all the 
desired properties in Section 4.2. If there are no unnecessary flat spots in 
r„ and if Xn is continuous or, more generally, if Disc(x„) C Disc(2;), then 
we can simply let fn = but more generally we cannot, and must insert 
extra flat spots in r„; we will return to that later. Since we will not pay 
close attention to the jumps, we can let fn be only a minor perturbation of 
r, obtained by inserting all necessary flat pieces, but only very short ones, 
so that the sum of all the lengths of the differences of f„ from r due to the 
addition of these flat pieces is less than where (5„ — ?• as n — )• oo. In that 
way, we can ensure that ||f^|| < ST and ||f^ — r'\\i — )• as n — )■ oo. 

We now construct the desired parametric representation {un,fn), again 
focusing on a single subinterval [s2^(e2)5 S2j+i(^2)]- We consider the case 
where it is necessary to add new flat spots to f„. Hence we will add extra 
flat spots to r in order to create f.„ over the subinterval. The endpoints of 
the interval correspond to points of increase in r by (4.11), so that we can 
let "Tn = r there. We will change r only in the interior of the subinterval 
[s+(e2),S2"j+i(e2)]. 

Suppose that we are considering inserting a flat spot in f„ at t. (We need 
to be careful because x and Xn could have countably many discontinuities 
in this region.) Let si be such that r(si) = t. Now we choose a second point 
S2 so that S2 > si) where, first, S2 is in the middle of a flat spot of r, so 
that r{s2) = t2 > t and, second, 1^2 — t\ and |s2 — si| are both very small. 
We do a new construction here. We remove a small portion / of the flat 
spot at level t2 and insert a flat spot of that length / at level t. In order 
to be consistent with the rest of constructed so far, to carry out this 
construction we move the entire function f„ to the right by the amount / 
between si and 52- Since we have deleted the interval of length / from the 
flat spot at t2, this step of the construction leaves the original function f„ 
completely unchanged outside the interval [si,S2]. Thus the change can be 
kept arbitrarily small. 

By this construction, we can make the Lebesgue measure of the set on 
which r is changed in the construction of f„ be less than any (5„ > 0, we can 
then let 5n — 0. At the same time, we can keep ||f^|| < 3T for all n > n*{e3). 
Hence we can achieve the desired Li convergence ||f^ — r'||i — )■ as well as 
II — '"II 0. In particular, it is easy to achieve over all m+ 1 subintervals 

(4.28) \\rn-r'\\i<ei forn>n*(e3). 
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We now consider the spatial portion of the parametric representation, Un 
where again we use the construction in Section 4.4. Over this subinterval 
[s2j{£2) , S2j_^_i{£2)], by this construction and the triangle inequality, we have 

\\un - u\\ 

— 1 1 ^n. ^ 1 1 

(4.29) + \\xn o f„ — X o f„|| + ||x o f„ — a; o r|| + ||x o r — u|| 

^ Jma,xi,Xfi) ~\~ \\Xji x\\ + Ws{x, Wf-n ''ID ~l~ 2JjYiax(^)) 

where Jmax and Wg are defined in (4.1) and (4.2), respectively, restricting to 
this subinterval [s2j(e2)) •s^_|_i(£2)]- 

In order to apply the bound in (4.29) to control the distance ||n.„ — n|| over 
[s2j{£2) , S2j_^_i{£2)], we let £4 be defined such that £4 < £3 for £3 in (4.23), 
and 

Ws{x,e4) < £1. 

We then let n** = n**{e4) = n**(£4,u,r, {xn}) be such that all of the follow- 
ing hold for n > n**: 

\\rn - r\\ < £4, 

(4.30) Jmax(2;„) < Jmax(2;) +£1, 

ll^n ^11 ^ 3t/jnax(2^) ~l~ £■!) 

where the uniform norms and the maximum-jump functions of x„ and x 
are restricted to the interval [527 (^2)5 S2j+i(^2)] and Jma.x{x) < £1 over this 
interval by the construction in Section 4.3. The last two relations follow from 
Lemmas 4.1 and 4.2. Combining (4.29) and (4.30), we have 

(4.31) ||u„-n|| <9£i <£ for all n>n**(£4). 

This final bound is the "weak link" in the collection of bounds we obtain for 
the three cases. Overall, if n > ri**(£4), then we obtain 

\\Un — u\\ V \\rn — r|| < £ 

with ||f[j|| < 3T. Combining (4.24), (4.26) and (4.28), we also obtain (over 
all of [0,1]) 

ll^^n - '"'111 < £2 +£1 +£1 < 3£l < £ 

for n > n**(£4). 
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